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ABSTRACT 

With the elementa ry energy release events introduced in a previous paper 



(iLiu &: Fletcher II2009I ) we model the chromospheric evaporation in flaring loops. 
The thick-target hard X-ray (HXR) emission produced by electrons escaping from 
the acceleration region dominates the impulsive phase and the thin-target emis- 
sion from the acceleration region dominates the low-energy thermal component 
in the gradual phase, as observed in early impulsive flares. Quantitative details 
depend on properties of the thermal background, which leads to variations in the 
correlation between HXR flux and spectral index. For lower temperature and/or 
higher density of the background electrons, the HXRs both rise and decay more 
quickly with a plateau near the peak. The plateau is less prominent at higher en- 
ergies. Given the complexity of transport of mass, momentum, and energy along 
loops in the impulsive phase, we propose a strategy to apply this single- zone en- 
ergy release and electron acceleration model to observations of flares associated 
with single loops so that the energy release, electron acceleration, and evapora- 
tion processes may be studied quantitatively. 

Subject headings: Sun: flares — Sun: particle emission — acceleration of particles 
— plasmas — turbulence 



1. INTRODUCTION 

The presence and nature of accelerated particles in solar flares are deduced primar- 
ily from secondary impulsive radiations. Because of the short thermalization and radiative 
timescales in the lower solar atmosphere, low-energy impulsive emissions from the dense 
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chromosphere to photosphere can be thermal as well. Gradual emission components are pre- 
sumably thermal and early in the flare are strongly linked to the chromospheric evaporation. 
The elec tron acceleration t i mescale is much sho rter than that for the chromospheric evap- 
oration (IMiller et al. 1119971 : lAschwanden II2002I ). and it appears that the two processes are 
decoupled. Ho wever, electron accele ration appears to be an inevitable component of the flare 
energy release (IKrucker et al. II2008I ). which itself is multi-scale as evident from the rich tem- 
poral, spatial, and energetic characteristics of flares, and electron acceleration in flares can be 
very efficient. Solar flares are caused by energy release on macroscopic scales, and the amount 
of nonthermal ele ctrons inferred frora observations can be a significant fraction of the back - 
ground electrons (IBrown fc Melrose 1 119771 : iFletcher &: Hudson I l2008l : iKrucker et al. 1 12009| ) . 
Energetic electrons in flares therefore are mostly accelerated from the ba ckground plasma. 



and the acceleration should depend in part on the background properties (jPetrosian fc Liu 



2004 : iLiu &: Fletcher I l2009l : iLiu et al. I l2009l ). One therefore may probe this dependence 



with observations over a relatively long timescale, when the properties of the background 
plasm a have changed signiflcantly. This is especially true f o r flares dominated b y single 



loops fIVeronig et al. I l2005l : iLiu et al. I l2006l : IXu et al. I l2008l : iHannah et al. 1 120081 ) , where 



the modiflcation of the background plasma through chromospheric evaporation is expected 
to affect the characteristics of nonthermal radiations. 

In the context of stochastic particle acceleration by turbulent electromagnetic fields, 
low-energy particles mostly couple with each other through Coulomb collisions. In the high- 
energy range the energy gain of particles through interactions with the turbulent fields 
dominates over the eri ergy loss through Coulomb collisions with lower energy particles and 
radiation (IBenz 1119771 ). A power-law high-energy particle distribution forms as this energy 
gain competes with the particle escape from the acceleration region through spatial diffusion 
with the power-law spectral index determined by the acceleration and escape rates. The 
particle distribution in the acceleration region generally has a thermal form at low energies 
and a high-energy power-law tail with the transition from the thermal to nonthermal com- 
ponent located at the energy Et, where the particle energy gain rate is comparable to the 
energy loss rate. Over the past few decades, it has been demonstrated that this mechanism 



is compatible with many fiare observations (IHamilton fc Petrosian 1119921 : iMiller fc Roberts 



19951 : iPetrosian fc Liu II2004I ). However, there are very few testable predictions from de- 



tailed modeling. The same is true for other particle acceleration mechanisms ( 



1985: Miller et al. 1997: Tsuneta fc Naito 1998: Aschwanden 2002: Drake et al. 



iolman 



20061 ). 



Given the important roles that turbulent electromagnetic fields play in the energy release 
processes of solar fiares, in a previous paper we pointed out that direct modeling of the 
dynamics of turbulenc e may advance our studies of fiares and lead to testable predictions 
(ILiu &: Fletcher 1 12009| ). In the simplest scenario, where the turbulence is generated in a 
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closed loop with a characteristic intensity and length scale — an elementary energy release 
event — the emission characteristics of the flaring loop are relatively well determined. We 
showed that these elementary energy release events might explain the soft-hard-soft spectral 
evolutioi i of hard X-ray (HXR) pu lses observed frequently both in chromospheric and coronal 
sources ( iBattaglia fc Benz 1 120061 ) . In this paper, we have a summary of the elementary 
energy release events (Section [2j) and quantify the chromospheric evaporation and X-ray 
emissions driven by this energy release (Section [3]). The model predicts faster rise and decay 
of HXR fluxes with a plateau near the turbulence peak for higher values of the ratio of 
to the thermal energy at this peak. The correlation between the nonthermal HXR flux and 
spectral index also depends on the temperature and/or density of the background electrons 
(Section H]). This simple model does not quantitatively address the plasma heating and many 
important transport processes in flaring loops. We show, however, that it is possible to use 
RHESSI observations to test the essence of the model and quantify the energy release, particle 
acceleration, and evaporation processes (Section [5]). Conclusions are drawn in Section [61 



2. Elementary Energy Release Events 



Elementary flare bursts were noticed in early spectral observatio ns of HXR flares, where 



these bursts are well correla ted among all available energy channels (jVan Beek et al. 1 1 19 74 



De Jager &: De Jonge IQTsI). Later high- r esolution observations reveal soft-hard-soft spec- 
tral evolution, which iGrigis fc Benz I (120041 ) attribute to elementary acceleration events. Al- 
though the duration and peak flux may vary signiflcantly from burst to burst, these events 
can be readily identifled by the soft-hard-soft behavior]^ Theoretically we called them el- 
ementary energy release events. The model for t he elementary energy release event and 



stochastic acceleration of electrons is described in iLiu fc Fletcher I (120091 ). For the sake of 



completeness, we summarize the key results here. The energy release is characterized with a 
turbulence generation length scale le and turbulence energy density 6g-B^/47r, where B is the 
background magnetic fleld and bo is a dimensionless parameter characterizing the fraction of 
magnetic energy dissipated during the event. Considering the energy equipartition between 
the turbulence kinetic and magnetic energies, the eddy speed at le is given by Veo = boVA, 
where va = -B/(47rp)^/^ is the Alfven speed and p is the background mass density. The 
characteristic timescale is then given by Tr = h/veo- 



^Decomposition of HXR light curves into elementary bursts i s not always possible, especially for 'ex- 
tended bursts', whose spectral evolution ca n be quite comphcated ( Hoyng et al. 19761 : Krucker et al. 20081 : 
Grigis fc BenTlbondlshao fc Huang Il2009[) . 
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Detailed analysis of c orrelation between HXR flux an d spectral index reveals compli- 

2008) ) ■ The shapes of these bursts also 



cated quantitative results (iGrigis &: Benz 



show significant variations (ISui et al. 1120071 : 



2004 



2006 



Li et al. 



20091 ) ■ These variations can be caused 



by changes in properties of the background plasma as the energy release triggers responses 
in coronal and chromospheric plasmas. The topology of magnetic fields may introduce addi- 
tional features. We focus on fiares associated with single loops so that particle acceleration, 
plasma heating, and chromospheric evaporation are the dominant physical processes in the 
impulsive phase. 



In the ID model for turbulence cascade and stochastic particle acceleration by lMiller fc Roberts 



( 119951 ) ■ the particle acceleration and plasma heating start almost instantaneously upon the 
arrival of turbulence at certain small scales, where background particles resonate effectively 
with turbulent plasma waves. Applying such a model to elementary energy events would 
imply almost universal rise profile of HXR bursts. There is little observational evidence 
for such a universal HXR rise profile, and the turbulence cascade and particle acceleration 
are truly 2D proce sses, which are highly nonlinear, complicated, and not well understood 
j.Tiang et al.lbood ). Here we use bB with b < to characterize the effective intensity of 
turbulent magnetic fields in resonance with background particles and assume that 



mm 



boVA/h 
= -b'^Vx/le 



for the rise phase , 
for the decay phase 



1^ 



where the 'dot' indicates a derivative with respect to ti me t and the Kraichnan phenomenol- 
ogy h as been used to obtain the turbulence decay time ( Kraichnan Ill965 : Miller &: Roberts 
19951 ). Both the kinetic and magnetic energy density of this turbulence is given by b'^pv\/2. 
The energy dissipation rate is then given by 



S(t) = 2pvib'/Td = 2pvib%A/l 



(2) 



Note that pv\ = B'^/Att does not change with time and the derivative of 6^ is 2bb. There are 
also uncertainties of factors of a few in these timescales. 

As in the standard quasi-linear theory, the wave-particle interaction rates are propor- 
tional to the effective turbulence energy density b'^pv\. In cases where the particle distribu- 
tion is isotropic, the pitch-angle averaged electron acceleration and scattering timescales by 
the turbulence can then be parameterized, respectively, with 



Sanb' 



SsTtb 



-2 



(3) 



where Tt = h/v, v is the electron speed, and Sa and Sg are two dimensionless coefficients 
characterizing the acceleration and scattering of electrons by the turbulent fields. We use 
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a leaky-box loss term to account for the spatial diffusion of electrons along magnetic field 
lines. The corresponding electron escape time from the acceleration region 



(4) 



where /n > L is the length of the energy release site. This equation is slightly different from 
that in lLiu fc Fletcher I (120091 ). Considering that the escape time cannot be shorter than the 
mean transit time lo/{y/2v), we include the second term on the right-hand side. This term 
is important only in the limit of weak scat tering when the scattering time Tsc is much longer 
than ~ /o/f (jPetrosian fc Donaghy Ill999l ). In a more accurate trea tment, one also needs t o 
consider the reduction of the scattering time by Coulomb collisions (jPetrosian fc Liu 1120041 ). 
The trapping of e lectrons by magnetic loo ps converging at both footpoints will also increase 
the escape time ( iMiller &: Roberts Ill995l ). In the regime, where scattering by turbulence 
dominates, the effects of Coulomb scattering, the tr ansit-time term, and ra agnetic trapping 
may be negligible and one obtains the expression in iLiu fc Fletcher I ( 120091 ). 



In general, the electron distribution function at a given time t depends on both the 
electron momentum and spatial location. For isotropic distributions, the spatially integrated 
electron distribution in the acceleration region only depends on the electron kinetic energy 
E: N{E,t) = No{t)g{E,t), where A^o = AIqu is the total number of free electrons in the 
acceleration region , n oc p is the average electron density, and A is the cross section of the 
flaring loop. We first ignore the structure of the flaring loop and assume that the evolution 
timescale of the distribution function is much shorter than r^. Then the normalized quasi- 
steady distribution function g{E,t) can be heuristically prescribed as: 



9iE,t) 



go{Et/kBT,,6) { E'/'exp{ 



1/2 



exp( 



-E/kBT,) 
-Et/kBT,){E/Et] 



for E <Et, 
for E>Et, 



where 



5 = 1/2 + (1 + r,,/4w)i/2 = 1/2 + {1 + SaSJl/[2'^'b%ikS, + 2'/%%)]}'/' , 
^o(r, 5) = [r^/^ exp(-r)/(5 - 1) - r^/^ exp(-r) + TT^^hri{r^/^) /2]-^ , 



(5) 

(6) 
(7) 



go is obtained from J g{E)dE = 1 and the time dependence is carried through Te(t), Et{t), 
and 6{t). § The transition energy between the thermal and power-law components 



Et = {27clnASJenY^^eyb 



^ A more accurate trea tment of g{E) with the effects of Coulomb coUisions taken into account can be found 
in lGallowav et al. I (|2005l ). and considering the loop structure will lead to a multi-thermal distribution, which 
may mimic a power law. 
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is obtained from the equality of E/Tac{E) to the Coulomb energy loss rate 

E,{E, t) = 23/27rlnAe%(t) /{Emy/^ (9) 

at Et, where Te, m, e, k-B, and InA ^ 20 are the electron temperature, mass, charge, Boltz- 
mann constant, and Coulomb logarithm, respectively: 

2^/^E'f%y{Sahm'/^) = 2'/WnAe'n/{Etmy/\ (10) 

The dependence of Et on b and n comes from the acceleration time and Coulomb energy loss 
rate, respectively. Then the (thermally) normalized transition energy 

r{t) ^ Et{t)/kBn{t) = ro[n{t)/no]'/'[bo/b{t)][no/TM (11) 

where, and in what follows, the subscript '0' indicates quantities at the peak of b{t). The 
dependence of r{t) on Tg is due to energy normalization to the thermal energy k-QTe. 

For convergence of the energy flux carried away by escaping electrons, the index of 
g{E)/Tesc, S — 1/2, must be greater than 2 (i.e., 6 > 5/2) if there is no high-energy cutoff in 
g{E). From equations ([nD, dSD, and (jl]), we then have 

r,, > 12re,, , b^ < [(1 + 2Sj3S,y/^ - l]kSJ{2'/Ho) = . (12) 

A is an upper limit for b^. As b^ approaches to A, the energy flux carried away by non- 
thermal electrons increases quickly near the turbulence peak. The damping of turbulence by 
nonthermal electrons will be important and equation ([1]) needs to be modified accordingly. 
These complex nonlinear effects ca n be considered wherever t here are sufficient theoretical 



and/or observational justifications (IBykov &: Fleishman II2009I ). There is also a physical up- 



per hmit on A. Since the turbulence energy density blpv\ should be less than the magnetic 
field energy density B^/Sti, we have 6o < A < 2~^/^. For bo > 2~^/^, the fiaring loop will not 
be well defined and the turbulence pressure may drive the acceleration site into an expansion. 
The model needs to be modified properly to account for these kinds of events. Therefore 
6o < A < 2^^/^ defines a self-consistent domain of the model. 

The escape of thermal electrons from the acceleration region can be treated as con- 
duction. However, scattering by turbulence will suppress the spatial diffusion of electrons 



near E'i, which will lead to a conductivity lower than the Spitzer conductivity (jjiang et al 



200' 



The loop structure a nd fiow motions may also affect the conduction significantly 



(lAntiochos fc Sturrock I Il978l ). In this paper, we do not seek more accurate treatment of 
electron transport at low energies. The formulae above is assumed to be applicable in the 
whole energy range. 
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From equations ([6]) and f|T2l) . we then have 



{(l + 2SJ[i5-0.5)'-l)S,]'/'-l}kSs 



where 



,2.c ^ {1 + 2./[(^- 0.5)2 -1]}^/^-! 
^ ^^'^^ = (l + 2./3)V2_i < 1 



One can show that ^(5, s) < 1 is equivalent to 5 > 5/2. It follows 



Sal Ss 



64A-4[i + SJ?,Ss - (1 + 25„/35,)V2] + 62A-2[(i + 25„/35,)V2 _ i] 

For Sa > 5„ e'(5,s) = {3/[(5-0.5)i^l]£^l5 = 1/2 + (1 + 3AV64)i/2, and 
{SaSs/l2f'^{le/lo) in agreement with llhi fc Fletcher I (|2o3). 



(13) 



(14) 



(15) 

1/2 

(16) 



3. X-rays from a Flaring Loop and Chromospheric Evaporation 

If we assume thin-target emission from the acceleration region (coronal sources) and 
thick target for escaping electrons (chromospheric footpoints), the overall specific intensity 
of photons observed at Earth is 

n^. - ^57w f {«(^- + /; £g;|d^ j ■ 

(17) 

where R = I AU, e, and o"(-E,£:) are, respectively, the distance between the Sun and 
Earth, the photon ener gy, and the angle-integrated bremsstrahlung emission cross section 



(IBethe fc Heitler Ill934l ). The first and second terms on the right hand side correspond to 
the thin-target acceleration region and thick-target footpoints, respectively, and transport 
effects from the acceleration region to the thick-target regi on have been ignored. ErjE^ t) and 
TesciE,t) are given by equations and (jl]), respectively (jPetrosian fc Donaghy Ill999l ). One 



can show that 4(^,t)w(^',t) = Ec{Eut)T,sc{Et,t)Et{t) / {EE'f^ = E^^{t)T,sc{t) /[{EE'f'^Tadt)] 
where we have used equation (fTOl) . From equation ([6]), we have Tescit) / Tacit) = 4/{[5(t) — 
l/2]2 - 1}. It follows that 

F(FtW (F' f\- 4[r(t)fcBTe(t)]2 

E,(E, t)rUE , t) - |[5(,)_i/2]2_i}(^^,)i/2 ■ (18) 



- 8 - 



Therefore I{e,t) only depends on N{E,t) and n{t), and the relation between the thin- and 
thick-target component is determined by N{E,t) alone, which itself is determined by r{t), 
n{t), and T{t). 

For elementary energy release events, equation ([1]) governs the evolution of b{t), which 
can be used to derive Et(t) and 6(t). One also needs to know the evolution of No(t), Te(t), 
and n{t) to obtain the evolution of N{E,t), which dictates the radiative characteristics. For 
events associated with closed loops, one may assume that Iq and A do not change with time. 
iVo(t) is then proportional to the average density n{t). In principle, one can solve the mass 



loops to st 


udy the detailed 


Klimchuk 


2006; 


Liu et al. 



20091 ). Here we are mostly interested in the spatially integrated properties. So a much 



simpler treatment is possible. 

One inevitable consequence of the impulsive energy deposition into closed coronal loops 
is chromo spheric evaporation. This e vaporation can be caused t hrough therrn al conduc- 
Antiochos fc Sturrock Ill978l ). non thermal particles (e.g . jBrown Ill973r ). and wave 



tion (e.g. 

or turbulence di ssipation at the footpoints (iFletcher fc Hudson I l2008l : iBrown et al. 1 12009 



Haerendel II2009I ). The details of these processes a re likely to be co mplicated and may affect 
low-energy emission characteristics significantly (IKlimchuk 1 12006| ) . For studies of X-rays, 
these processes may introduce short timescale features, which are energetically less impor- 
tant. This paper focuses on the dominant X-ray emission component on a timescale of 
Tr = le/h^vx- For the sake of simplicity, we assume that the average density increase rate in 
the coronal loop is proportional to the energy dissipation rate S, i.e., 

h{t) = C,nit{t)/{fnjl) = 2CenibW/le , (19) 

where Cg is a dimensionless coefficient and the subscript indicates quantities at the flare 
onset as is in what follows. Note that pvj^ = B^/An gives the normalization of the energy 
density and does not change with time, rii can be considered as a model parameter to be 
determined by observations. 

The above energy conservation equation can be generalized as follows 

3kB[nit)f{t) + n(t)T(t)] = 3kBCeT{t)ni[t{t) - Ait)]/{pvl) , (20) 

where we have assumed that the electron and ions temperature are equal and the internal 
energy density is given by SnkBT with T{t) = Te{t). A(t) represents the overall energy loss 
rate from the flaring loop. If energies in forms other than the background magnetic fleld, 
turbulence, and internal energies of the thermal plasma are negligible, we have ?)kB[n{t)T{t) + 
n{t)T{t)] = t{t) - A(t), which gives 

Ce = pvl/{3kBTni) . (21) 
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Therefore CgT should be approximately a constant and we recover equation f|T9l) for T = 
and A = 0. 



Equations ([T]) and f[T^ then give 



1 + Cem^ - bt]/{2bo) ^a^ + C,by{2bo) 
a^ + Cyj2 + C,[bl-b\t)] 



for the rise phase , 
for the decay phase . 



(22) 



Combining equations fl2^ and ([T]), we have for the rise phase 

I *f/\(^; I U;|/fc7/ /till ~ 



.■1/21 



h(t) + alnJrioW^ In <^ ^^^^ — 

2(6^ - bl)/bl = hit) 

hit) - 2a,[n,/no]'/Htg-'[n/in,al)Y/' - tg-'[no/in,al)Y/'} 



for 1 > Cebf/i2bo) 

for 1 = Cebt/i2bo) 
for 1 < Cebf/i2bo) 
(23) 

where Oi = [Ce6f/(26o) - t^Ao = 5/(47^po)^/^ and hit) = 2[nit)^/^ - nj^]/nj\ From 
the second expression, we see that the flare starts after tg = — /e/(2fAo^o) when bits) = 
and nits) = 0. The third expression can also be extrapolated backward in time to a point 
ts when nits) = 0, i.e., Ce6f/(2feo) = 1 + Ceb^ i2bo), where bs = bits)- Note that n^/ni = 



+ Cebl/2 = C,bl/2 



2botsVAo/h = -1 



1 + Ce(&o bf)/i2bQ), and from the third expression we have 



Cebt/i2bo) 



nl/2 



l + C,ibl-bf)/i2bo) 



nl/2 



l + CM-bf)/i2bo 

= -1 + iibo/bs)' - ir'^hg-' [ibo/bs)' - 1] 

The solution for the decay phase can be obtained numerically. 



Cebt/i2bo) 

1/2 



(24) 



We consider first the solution in the middle of equation (!23|) . Then Cg = 2bo/bj = 
2no/inibl), and n = nob^/bQ and no(l + 2/bo — 26^/6q) for the rise and decay phase, respec- 
tively. If we assume a coronal background level for 6, and n,, the density increases by a factor 
of &o/^f 1 + 2/^0 in the rise and decay phase, respectively. There is more evaporation 
for stronger events especially during the rise phase. From equation (|T9l) . one has the energy 
dissipation in the rise and decay phase of (1 — bj/bQ)ib^/2)pv\ and blpv\, respectively. The 
energy dissipated in the decay phase is equal to the total turbulence energy at the peak as 
expected. In the rise phase, equation ([T]) shows that there is a continuous energy deposition 
into the loop and the energy dissipation depends on the evaporation process. 



4. Results 



The solution of equations ([T]) and fl22l) critically depends on the sign of 1 — Cebf/2bo. 
For Cebf = 2bo and given bo/ A, Sa/ Ss, tq, uq, and Te(t) = Tgo, NiE,t) is determined for 
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Fig. 1. — Evolution of g{E,t)n{t)/no during the rise (upper) and decay (lower) phase for 
To = 4 (left) and 8 (right). The time unit is tq = le/{boVAo)- Cgbf = 2bo, Sa/ Ss = 100 and 
A = 2^/^6o for both cases. 



the time in units of tq = le/{boVAo)- We will first study this simple case with remaining 
constant and discuss more reahstic applications of the model to flare observations in the next 
section. For bo = 2~^/^A, and Sa/Sg = 100, Figure [H shows the evolution of g{E,t)n(t)/no. 
It is evident that, while the electron distribution is thermal and nonthermal, respectively, at 
low {E <C ro/cB^'e) and high {E ^ ro/cB^'e) energies, in the intermediate energy range {E ~ 
ro/cB^e) it experiences a transition from nonthermal to thermal. Moreover, the evolution 
of the gradual thermal component is similar for ro = 4 and 8, in contrast to their distinct 
impulsive nonthermal components. 

Figure [2] shows the time dependence of g{E, t)n{t) / hq at several energies. In combination 
with Figure [H it is evident that, for the ratio of transition energy to thermal energy at the 
turbulence peak ro = 4, electrons at 6 ksTe are nonthermal and impulsive near the peak of b 
and become thermal in the gradual decay phase. At = 12 and 25 k^Te, the distribution is 
nonthermal and impulsive. For Tq = 8, electrons at 6 k^T^. are nonthermal in the early rise 
phase and become part of a thermal distribution be fore the peak o f b. Its sharp rise near 
the peak of b mimics the well-known Neupert effect (iNeupert 1 119681 ) . Electrons at 12 /cB^^e 
become thermal in the decay phase. Compared with the model with tq = 4, the most distinct 
feature is the much sharper decrease of g at E = 12 and 25 /cb^e in the early decay phase. 
This is mostly caused by the relatively high values of Et/k-^Tf., which makes the density of 
nonthermal electrons very sensitive to Et and therefore b{t) and n{t). 
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4 

i -1 -6 -5 -4 -3 -2 --8 -7 -6 -5 -4 -3 -2 

Logio[Eg(E)n(t)/n(0)] Logio[Eg(E)n(t)/n(0)] 

Fig. 2. — Evolution of g{E,t)n{t)/no (upper) and correlation between 6 and g{E,t)n{t)/no 
at E = 25fcB^e (lower) of the evaporation models shown in Figure [H The dashed and dotted 
lines in the lower panels correspond to the rise and decay phases, respectively. 



Figure [3] shows the corresponding evolution of el{e, t). To increase the dynamical range 
of the plot, we show the energy spectrum el{e, t) instead of the photon number flux spectrum 
I{e,t)/e normally plotted in papers on specific observations. The thin and thick lines are 
for the first and second term of equation ( fT7|) . and the upper and lower panels are for the 
rise and decay phases, respectively. The distinctions between the thermal and nonthermal 
components are less prominent than those for electrons primarily due to the dominance of the 
thick-target emission and contributions of high-energy electrons to low-energy X-rays. The 
nonthermal HXRs are always dominated by the thick-target component, which dominates 
the rise phase of the tq = 4 model as well. The thin-target component only dominates the 
low-energy thermal emission in the late rise and decay phases. It is interesting to note that 
a low-energy spectral flattening in the total spectrum can result from competition between 
the thin- and thick-target components, whi ch provides an alternative to th e return current 
explanation for this observed phenomenon (jZharkova fc Gordovskyy II2006I ). 

The upper panels of FigureHlshow light curves of el{e, t) at several energies. The impul- 
sive phase is mostly dominated by the thick-target emission and the thin-target component 
dominates in the gradual phase especially at low energies. Compared to the case with tq = 4, 
the HXR fluxes of rg = 8 both rise and decay more sharply with a plateau near the peak. 
This plateau is less distinct at higher energies. The dependence of HXRs on tq is a direct 
consequence of the dependence of the nonthermal electron density on tq. A high value of 
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"5.0 0.5 1.0 1.5 2.0 §.0 0.5 1.0 1.5 2.0 



LogioCeAsr,,) LogioCeAsT'e) 



Fig. 3. — Evolution of the energy spectrum el during the rise (upper) and decay (lower) 
phase for tq = 4 (left) and 8 (right), where I is the specific intensity rather than the photon 
flux, as is more often used. The unit is 2^/^ae'^/o^'^o(^B^e)^''^/(37rm^/^c^i?^), where a is the 
fine structure constant and c the speed of light. The thin and thick lines are for the thin- 
and thick-targets, respectively. 



ro can be caused by a low temperature and/or high density of the background electrons for 
a given acceleration timescale. The model therefore predicts that sharper rise and decay of 
HXR emission imply lower temperature and/or higher density of the thermal component. 
This is the key prediction of the model. 

Since the evolution of 5{t) is determined by h{t) (eq. [B]), the dependence of the HXR 
light curves on Tg and n will affect the correlation between HXR flux a nd differential photon 
spectral index 7 = 1 — d l nJ/dlng frequently observed for HXR pulses (iGrigis &: Benz 112004 



Battadia fc Benz Ibood : ICriris fc Benzlbood . boosh . (Note that since 7 is usually defined 
in terms of photon number flux spectrum instead of energy flux spectrum, it differs from the 
index of / by 1.) The lower panels of Figure H] show how the correlation changes from the rise 
to decay phase at e = 25fcBT'e- The corresponding correlation between 5 and g{E ,t)n{t) /tiq 
is shown in the lower panels of Figure [2l Since the thick-target emission does not depend 
on the density of the target region, the correlation between 5 and g{E,t)n{t) is similar to 
the correlation between 7 and el. The thin-target emission scales with n{tYq(E, t), and 
el therefore decays more slowly than n{t)g{E,t). iGrigis fc Benz I (120041 ) claims that the 
correlation for these pulses "follows tendentially a slanted V, with the rise phase forming 
the flatter leg". This is in agreement with the ro = 4 case. For higher values of ro, our 



- 13 - 



o 

to 
o 

J - 



bO -6 

o 

hJ -0.: 



6 



b/A 
6k,T, 
12 is TV 

25 kaT, 



t/To 



Logio(el) 



o 

< - 



o 



b/A 

llk,T, 
25 tj J' 



t/To 



-4 -2 



Logio(el) 



Fig. 4. — Light curves (upper) and correlation between differential photon ffux index 7 and 
eJ at e = 2bk-BTe (lower) of the evaporation models shown in Figure [3l The dashed and 
dotted lines in the lower panels correspond to the rise and decay phases, respectively. And 
the thin and thick lines are for the thin- and thick-targets, respectively. 



model also predi cts an opposite evolu tion with the decay phase forming the flatter leg. Such 
pulses do exis t (IGrigis fc Benz 1120081) . The strongest pulse of the 25-Feb-2002 flare shown 
in Figure 4 of iGrigis fc Benz I (120041 ) is one example. It is also interesting to note that this 
pulse has a relatively gradual rise and sharp decay phases in agreement with our model 
prediction at high ener gies. The majority of pulses from the 27-June-1980 flare studied by 



Lin &: Schwartz I (119871 ) also seems to be compatible with the latter model as shown in their 
Figure 8. 



5. Discussion 

We have shown that the model has many characteristics similar to t he observed ther - 



mal and nonthermal X-ray emissions, especially for early impulsive flares (ISui et al. 1120071 ). 
It also predicts that variations in properties of the thermal background plasma result in 
significant differences in the nonthermal particles. To test the model more quantitatively 
with observations, one needs to address the heating of electrons and ions and other impor- 
tant processes in flaring loops. The evolution of Tg and ion temperature Tj depends on the 
ill-understood plasma heating by turbulence and chromospheric evaporation. The inhomo- 
geneity along the loops and related production of plasma waves will also affect the thermal 
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energy gain of background particles (IZharkova fc Gordovskyy II2006I : iKontar fc Reid II2009I ). 
In cases where the energy release site is localized at the looptop regions, one also needs 
to consider transport effects from the acceleration region to the thick-targ et region. Some 



transport processes may modify the electron distribution by inducing waves (IKontar fc Reid 



20091 ). which may propagate into the acceleration region and affect both the acceleration and 
transport processes. The trapping of energetic electrons in the accel eration region also mod- 
ifies the resulting photon spectral evolution (IGrigis fc Benz 1 12006| ). The problem is even 
more complex when one consid ers the effect of c onduction and radiation from coronal loops 
and chromospheric footpoints (iLiu et al. II2009I ). For strong evapora tion, shocks ma y form 
making the large-scale bulk motion an important energy component (ILi et al. II2009I ). 



There are already extensive numerical studies of the d ynamics of flaring loops (ILi et al 



19971 : lAUred et al. II2005I : iKlimchuk II2006I : iLiu et al. II2009I ) . Significant uncertainties remain. 
The essential difficulties are rooted in the highly dynamical nature of turbulent plasma in 
a strongly magnetized environr nent and we still lack effective tools to address this phe- 
nomenon in quantitative detail (j Jiang et al. 1 120091 ). It is therefore not surprising that flare 
HXRs often show features over a broad scale range. For some flares, features from the du- 
ration of t he observed HXRs down to the time resolution of the instrument can be readily 
identified (lAschwanden et al. Ill998l ). Even for the relatively simple early impulsive flares 



with the HXRs dominated by a few distinct bursts (ISui et al. 1120071 : iLi et al. II2009I ). finer 



structures with a tiraescal e much shorter than the duration of individual bursts are expected 
(lAschwanden et al. Ill996l ). Flare HXRs are therefore an intrinsically complex phenomenon. 
Detailed deductive modeling of individual flares is deemed to be tedious if possible, especially 
for extended bursts without distinct dominant scales. 

On the other hand, for simple flares with the HXRs dominated by a few distinct bursts, 
which usually follow the soft-hard-soft evolution, our model can be readily applied to in- 
dividual bursts to account for these energetically dominant components. The energetically 
less important residuals from a successful model fitting to specific flare observations can be 
attributed to small-scale processes ignored (or averaged out) in the model. An expansion 
approach in energetics therefore may be taken to achieve more and more quantitative mod- 
eling. The elementary energy events proposed here therefore may reflect more about the 
fact that some flares are relatively simple and may be modeled in quantitative details rather 
than some fundamental physical processes. If there are universal processes operating in all 
flares, these simple flares will be good candidates for us to disentangle these processes by 
modeling in detail, which may lead to interesting testable predictions. With the elementary 
energy release events proposed here, one can carry out detailed modeling of early impulsive 
events, where properties of the background electrons can be readily obtained from RHESSI's 
high resolution observations of X-ray emissions. Quantitative tests of the model are possible. 
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especially for the study of energy release during the impulsive phase when one may assume 

Ti = Tf. = TeQ. 

The energy loss needs to be treated properly to apply the model to observations espe- 
cially for the gradual decay phase. To partially take into account the effect of energy loss 
during the impulsive phase, one may replace equation (fT9l) with 



n{t)=C,n,[t{t)-A]/{pvl) 



(25) 



where A = KrU^ /tI^"^ is the radiativ e cooling rate and /t^ — 1.42 x 10^^^ ergs cm^ s"^ K^/^ 
(iTandberg-Hanssen &: Emslie Ill988l ). Figure [5] shows how the density variation is reduced 
due to inclusion of the radiative cooling for Aq/Sq = 0.22/6q, Ceiiibl = 2nQ, and bo = 1. 
The density decrease for the result with radiative cooling considered may not represent 
observations accurately because depletion of plasmas from coronal loops produced by flares 
usually occurs when the plasma temperatu re is already so low that there are no significant X- 
ray emissions (jWarren fc Antiochos 1120041 ). The temperature evolution must be considered 
before the density reaches the peak value. This can be done when applying the model to 
individual flare observations. 



without radiative cooling 
with radiative cooling 




Fig. 5. — The evolution of the average density n with (solid) and without (dashed) the 
radiative cooling taken into account. See text for details. 



Flares associated with single loops have been studied extensively by many author s 
JVeronig et al. Il200,4 ILui et al. Iboonl : ISui et al. \\m)^ . Iflannah et allboosl : IXu et al. hoO^ . 
Two flares occurring on 13-Aug-2002 after 5:30 UT are also good examples ( jLi et al. II2009I ). 
Among these fla res, there are events with the HXR emission dominated by a single pulse 
(ISui et al. 1120071 ). Our elementary energy release events can be applied to these flares di- 
rectly to study all aspects of the impulsive phase. These pulses usually have a duration of 
less than 1 minute. It may be difficult to obtain the detailed spectral evolution, however. 
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one can use the model to reproduce the observed hght curves with a prescribed temperature 
evolution. Te, Ce, and uq determine the thermal X-ray light curves. The photon spectral 
index at the HXR peak can be used to derive bo/ A, which determines the evolution of the 
HXR spectral ind ex. The time unit tq can be adjusted to fit the duration of the HXR pulse 



(ILee et al. 1 119951 ). The shape of the HXR pulse can be fitted by adjusting Sa/Sg and tq. 
The overall normalization Nq depends on the size and density of the looptop acceleration 
region. Therefore a good fit to such flares would lead to well-defined model parameters. If 
observations can also give a reasonable measurement of the source size Iq, from equations ([H]) 
and f|T2l) . one can show that the model parameters uq, Tgo, bo/ A, Sa/ Ss, tq, Iq must satisfy 
the following equation 

ro^BTeo = {2'/WnAnoloSa/Ssy/'e'{A/bo)/[{l + 2Sj3S,y/' - if' . (26) 

l^Sg cLiid l^Sfi are introduced in equation ([3]) and can be obtained in term of A with equations 
(fT2|) and (Q. The magnetic field B needs to be determined from other means to get le (also 
in term of A) from tq and bo/ A. It is evident that A can not be derived from the above 
model parameters. By replacing A with 1^, one, however, can express Sa and Sg in terms 
of le, which are intimately related to the energy release and electron interactions with the 
turbulence and should not change significantly from flare to flare. Of course, if the turbulence 
speed &o^A can also be determined through other means, one can get le from tq and then 
derive bo, A, Sg, and Sa- The capability of the model to reproduce observed light curves with 
the constraint fl2^ can be considered as a test of its validity. From fitting of a large enough 
sample with different peak HXR spectral indexes, one may obtain statistical measurements 
of Ce, le, bo. A, Sa, and Sg. 

There are also flares with relatively long HXR pulses 1 min), e. g ., the 20-Sep-2002 



flare starting at 9:23 UT and the flare studied in detail by iLiu et al. I (120061 ). Given the 
relatively gradual HXR peak of these events, the elementary energy release events will not 
produce a good fit to the light curves. For these events, one can obtain detailed spectral 
evolution. Equation ([6]) therefore can be used with the observed 6{t) to derive the evolution 
of b{t) and S(t) during the HXR pulses. The temperature evolution Te{t) can also be 
obtained from the spectral evolution. With Te{t) and 6{t) as inputs, one can again fit 
the observed light curves with the model to test its validity, especially the evaporation 
model described with equation (fT9|) . In a more accurate treatment, one should replace 
equation (fT9!) with an appropriate energy conservation equation to have more insights of 
the evaporation a nd energy release proce sses. For more general cases, where there can be a 



preheating phase (IBattaglia et al. I l2009l ) and/or many HXR pulses (jVeronig et al. II2005I ). 
suggesting continuous energy release, one can treat bit) as an input and apply the model to 
observations to derive its evolution, which can be used for further study of plasma heating, 
particle acceleration, and chromospheric evaporation. 
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6. Conclusions 

Although the particle acceleration timescale is short, the energy release processes of 
solar flares are multi-scale phenomena and some long timescale processes may have imprints 
on the characteristics of nonthermal components. In this paper, we demonstrate how the 
evolution of X-ray emissions may be affected by the chromospheric evaporation in the context 
of stochastic acceleration by turbulent electromagnetic fields. We consider the simplest 
scenario, where the energy release in closed loops is described with a turbulence intensity 
boB and generation scale le- It is shown that the HXR flux has a sharper decay phase for 
events with lower temperature and/or higher density of the background electrons, which 
makes the transition energy between the thermal and nonthermal component much higher 
than ksTe. As the ratio, ro, of the transition energy to the thermal energy at the HXR 
peak is varied from low to high values (by varying the density and temperature) the HXR 
behaviour changes between two types of spectral evolution. If this ratio is low, the flux 
change in the rise phase is higher than in the decay phase for a given change in the spectral 
index (the commonly observed behaviour) and vice versa if Tq is high — a behaviour which 
is also observed, though less frequently. 

Solar flares involve many physical processes that give rise to signiflcant uncertainties in 
our understanding of the plasma heating and particle acceleration. Even for relatively sim- 
ple flares associated with single loops, the mass, momentum, and energy transport along the 
loops during the impulsive phase can be highly dynamical and complicated. It is therefore 
challenging to study the flare energy release with all the relevant processes treated properly. 
The problem can be simplifled significantly if we are mostly concerned with the energeti- 
cally dominant processes such as the energy release, particle acceleration, and chromospheric 
evaporation averaged over a relatively long timescale of observational interest. We discuss 
how the model may be applied to RHESSI observations of flares for quantitative investiga- 
tions, especially for early impulsive flares with the HXRs dominated by a single pulse with 
soft-hard-soft spectral evolution. 

While the short acceleration timescale may suggest microscopic processes, the large 
amount of energy released during flares reveals a predominantly macroscopic process. One 
critical aspect of solar flare studies is to understand the connection between the microscopic 
particle acceleration and the macroscopic energy release. The latter has been studied with 
magnetohydrodynamical simulations. The study of particle acceleration is limited to a few 
very speciflc mechanisms. Although the high degree of freedom on large scales (in connection 
with the magnetic field configuration and initial and boundary conditions) inevitably lead 
to multi-scale features giving rise to significant complexities, the small scale processes are 
expected to follow some patterns as least in a statistical sense. One therefore may focus on 
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flares with relatively simple HXR light curves. The model presented here makes it possible to 
infer the energy release rate from the observed HXR emission for these flares. It is therefore 
possible to probe the macroscopic energy release processes with the impulsive nonthermal 
emission. Successful applications of the model to flare observations will lead to quantitative 
measurements of the particle acceleration and scattering rate by turbulent electromagnetic 
fields. The obtained b{t) can be compared with MHD simulations of fiares or models of large 
scale energy release processes for more quantitative and self-consistent modeling. 

This work is supported by the EU's SOLAIRE Research and Training Network at the 
University of Glasgow (MTRN-CT-2006-035484), the Ministry of Science and Technology 
of China (Grant No.2006CB806302) and by RoUing Grant ST/F002637/1 from the UK's 
Science and Technology Facihties Council. 
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